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Using a cylindrical soliton solution to the four-dimensional vacuum Einstein equation, we study 
non-linear effects of gravitational waves such as Faraday rotation and time shift phenomenon. In the 
previous work, we analyzed the single-soliton solution constructed by the Pomeransky’s improved 
inverse scattering method. In this work, we construct a new two-soliton solution with complex 
conjugate poles, by which we can avoid light-cone singularities unavoidable in a single soliton case. 
In particular, we compute amplitudes of such non-linear gravitational waves and time-dependence 
of the polarizations. Furthermore, we consider the time shift phenomenon for soliton waves, which 
means that a wave packet can propagate at slower velocity than light. 

PACS numbers: 04.50.+h 04.70.Bw 


I. INTRODUCTION 

A lot of gravitational solitons in general relativity, which describe gravitational solitonic waves propagating through¬ 
out spacetime, have been found within the framework of so-called inverse scattering method Cli- In particular, it has 
attracted a lot of relativists since it can generate black hole solutions in an axisymmetric and stationary case d, 0 , 
in addition to exact solutions describing non-linear gravitational waves on various physical background. The method 
was immediately generalized to the higher dimensional Einstein equations (see the references in Q) but such a sim¬ 
ple generalization to higher dimensions tends to lead to singular solutions. However, 10 years ago, Pomeransky (Ef 
succeeded in modifying the original inverse scattering method [6| so that it can generate regular solutions even in 
higher-dimensions. Thanks to this, as for five dimensions, several black holes solutions were found Hi- 

A diagonal metric form of a cylindrically symmetric spacetime makes the vacuum Einstein equation to an extremely 
simple structure of a linear wave equation in a fiat background. As a solution of such a equation, Einstein-Rosen wave, 
which can be interpreted as superposition of cylindrical gravitational waves with a + mode only, were obtained [10- 
However, the existence of nontrivial non-diagonal components of a metric drastically changes the structure of the 
Einstein equation since it generally yields a x mode together with non-linearity. Piran et al. [H numerically studied 
non-linear interaction of cylindrical gravitational waves of both polarization modes and showed that a + mode converts 
to a X mode. This phenomenon is named gravitational Faraday effect after the Faraday effect in electrodynamics. 
Tomimatsu [9[ studied the gravitational Faraday rotation for cylindrical gravitational solitons generated by the inverse 
scattering technique. Moreover, the interaction of gravitational soliton waves with a cosmic string was also studied 
in |1Q|-[i^. As one of new attempts to understand strong gravitational effects we have recently constructed a new 
cylindrically symmetric single-soliton solution from Minkowski seed by the Pomeransky’s inverse scattering method 
and clarified the behavior of the new solution including the effect similar to the gravitational Faraday rotation. 


* tomizawasny@stf.teu.ac.jp 
1 tmishima@phys.ge.cst.nihon-u.ac.jp 



2 


In this paper, for further investigation we first construct more complicated gravitational two-soliton solutions with 
complex conjugate poles by the Pomeransky’s method, which are considered as cylindrically symmetric gravitational 
waves with rich structure and can be used to study gravitational non-linear effects. One of important and remarkable 
features of the solutions is that there are no null singularities which appear generally in the case of single-soliton 
solutions. We may therefore adopt the picture similar to the scattering theory that from the past null infinity one 
gravitational wave packet comes into ‘the interaction region’ near the symmetric axis and after reflection leave the 
region for the future null infinity. The behavior of the wave packet can be analyzed by following the time sequential 
images and also by comparing the physical quantities measured in the past and future infinity. As an interesting 
example, the change of polarization of two independent modes will be treated with the both methods 

In the next section, we present the Kompaneets-Jordan-Ehlers form 0 in the most general cylindrical symmetric 
spacetime and the useful quantities (amplitudes and polarization angles) for analysis of non-linear cylindrical gravita¬ 
tional waves, which were first introduced by Piran et a/.[ 8 |] and Tomimatsu [9[ . In Sec. lIIIi using the inverse scattering 
method improved by the Pomeransky [5| , we will generate a two-soliton solution with complex conjugate poles from 
Minkowski spacetime. In Sec JIVI we will analyze the obtained two-soliton solution by computing the amplitudes and 
polarization angles for ingoing and outgoing waves. In this section, in particular, by seeing time-dependence of po¬ 
larizations, we study gravitational Faraday effect. Furthermore, we will mention the difference from the single-soliton 
solution in our previous paper, and moreover clarify the difference from the Tomimatsu’s two solution [9[ generated 
by the Belinsky-Zakharov’s procedure. In Sec. [Vl we will devote ourselves to the summary and discussion on our 
results. 


II. FORMULAS 

We assume that a four-dimensional spacetime admits cylindrical symmetry, namely, that there are two commuting 
Killing vector fields, an axisymmetric Killing vector d/d(j) and a spatially translational Killing vector djdz^ where the 
polar angle coordinate 0 and the coordinate z have the ranges 0 < 0 < 27r and — oo < 2 ; < oo, respectively. Under 
these symmetry assumptions, the most general metric can be described by the Kompaneets-Jordan-Ehlers form: 

ds^ = e^'^{dz -h ujd(l)Y + p^e~^'^d(j)^ -\- — dt^), (1) 

where the functions uj and 7 depend on the time coordinate t and radial coordinate p only. Following Ref. [H, [9[, 
we introduce the amplitudes 
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where the advanced ingoing and outgoing null coordinates u and v are defined hy u = {t — p)/2 and v = (t -h p)/ 2 , 
respectively. The indices and x denote the quantities associated with the respective polarizations. Then, the 
vacuum Einstein equation can be written in terms of these quantities, actually, the non-linear differential equations 
for the functions ip and uj are replaced by 
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The ingoing and outgoing amplitudes are defined by, respectively, 

^ = ^A\ + Al, B = + (12) 

and the polarization angles Oa and Ob for the respective wave amplitudes are given by, 

A B 

tan26>^ = —tan26>B = (13) 

III. TWO-SOLITON SOLUTION 

In this work, as a seed, we consider Minkowski spacetime written in cylindrical coordinates whose 2x2 part of the 
metric go := {goab) (<^ 5 ^ = and metric function /o := are given by, respectively, 

go = diag (l, p^) , /o = 1. (14) 

Following the inverse scattering method which Pomeransky improved we construct a two-soliton solution with 
complex conjugate poles. First remove trivial solitons with (1,0) at t = ai and t = a 2 { 0.2 = cii) from the seed metric 
and then we obtain a metric 


g'o = diag 






(15) 


where w := pi/p {pi = JI 2 = \/{t — diY ~ ~ — ai)). Note that 0 i is a complex parameter and the bar denotes 

complex conjugation. Next add back non-trivial solitons with (l,a) and (l,d). Then we obtain two-soliton solution 
as 


^ ^ (flo)acTOc^^(r ^)kim^Ji\go)db 
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f — f det(rw) 
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5'o(p, t, A) is a generating matrix for the seed g^ which is given by 

(p 2 + 2t\ + A 2)2 


dip A A) = diag 


(/il-A) 2 (/l 2 -A)' 


, p2 -(- 2tA + •> 


(18) 

(19) 


( 20 ) 


where pk = P^/Pk {k = 1 , 2 ). 

We present the metric in the Kompaneets-Jordan-Ehlers form [^, which describes the most general cylindrically 
symmetric spacetime, as 

ds2 ^ ^ ^ ^2g-2V-^^2 ^ (21) 

where the functions t/;, uj and 7 depend on the time coordinate t and radial coordinate p only and they are explicitly 
given by 
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with 


A = 2^ 
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^ 9 ’ 

c = 2aai. 

Here, 5R[ ] denotes a real part of [ ]. 

After using the time translational invariance of the system to rewrite the parameter ai as iq (i is and g' is a 
positive number), we can simplify the metric by introducing the following new coordinates (x^y): 

t = qxy, p = q\/{x‘^ + l){y^ - 1). (31) 

Note that 


(28) 

(29) 

(30) 


Hi= Hi= q{x + i){l -y). 

Let us put a = ttr + ttii (a^, Uj € R). In the coordinate system {x,y), the metric can be written as 


(32) 


.2 y f, 2^ ,,2 ^ ( dx^ dy^ 

ds _ -i^dz + -d<l>j +p-d<l> + 4ogg^4(^2 + y2)5 + ^i; ’ 


(33) 


where the metric functions A, y and Z are 


= o-iiy - lf{y + 1)® + 2af(y + 1)^ {al{y - lf{y + 1)^ + 64^^ {x'^{y{9y - 8) + 1) + 2x‘^{y{y + 4) - 3)y^ + y® + y^)) 
-512aiOry^a;(y + 1)^ {x^ - (y - 2)y) (a;^(2y - 1) + y^) + aj{y - lf{y + 1)® 

+128alq\y + if (2a;® + xfiS - 3y)y - 1) + 2x‘^y\2{y - 2)y + 3) + y® - y^) + 4096y^ {x^ + y^ , 
y = + 2af {if - l) {al (y^ - l)^ + 64^^ (a;^ (9y^ “ l) + 2a;^ (y^ + 3) y^ + y® - y^)) 

-1024aiary^a; {x^ + l) y (y^ - l) (a; - y){x + y) + (y^ - l)"^ 

+128a^y^ (y^ - l) (2a;® + a;^ + (4a;^ + l) y^ - 3 (a;^ + 2) a;^y^ + y®) + 4096y"^ (a;^ + y^ , 

Z = -32q^(y + l)(afa;(y - l)(y + 1)® (a;^(l - 3y) + (y - 3)y^) + afor(y - l)(y + 1)® (a;"^ - 3a;^(y - l)y - y®) 

+aia;(y - 1) (al{y + 1)® (a;^(l - 3y) + (y - 3)y^) + 64y^ (a;^ + y^)^) 

+ar (al{y - l)(y + 1)® (x^ - 3x^(y - l)y - y®) - 64y^ (x^ + y) (x^ + y^)^) . 


IV. ANALYSIS 

Let US investigate how the new gravitational solitonic waves propagate throughout spacetime from the several 
viewpoints. For convenience of explanation, we introduce the modulus k and the angle 0 of the complex parameter a 
defined by eq. (33), 


k = |a|, 0 = Arg(a). 

In the following analysis, we only consider the case of g' = 1, because the parameter q can be normalized by the scaling 
of coordinates. 
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At the beginning we show the qualitative behavior of the waves near the axis briefly. The figures n = 0, • • •, 7 
in FIGil] display the various behaviors of total amplitude (Atot = +5^) from the spacetime viewpoint, which 

are arranged from left to right in each figure, respectively. The figures are plotted in the range of — 5 < t < 5 and 
0 < p < 5 under the parameter-setting g = 1 , p = 2 and 0 = nir/4 {n = 0 ~ 7). Some of the behaviors are also 
displayed in FIGl 2 ]by superimposing the instantaneous graphs that correspond to t = 0 and t = (n = 1 ^ 3), 
respectively. From the behaviors of the amplitudes we may consider the gravitational solitonic waves to be regular 
wave packets, which first come into the region near the symmetric axis from past null infinity, and leave the axis after 
reflection for the future null infinity. On first inspection the various complex behaviors seem to be generated by a 
certain kind of non-linear effect near the axis. For example, shown in the case of n = 1 in FIGHJand FIGl 2 l at the 
angle ^ = 7 r /4 the initial two wave packets coalesce into one packet, while at the angle 0 = 37 r /4 the initial single 
wave packet splits in two. These remarkable phenomena seem to occur very near the axis, so that we may say that 
the neighborhood of the axis is the region where non-linear effects become highly strong. 

In the subsequent subsections, we will see the detailed behavior of the waves propagating near the boundaries of 
spacetime, particularly focusing on its non-linear effect. 




FIG. 1: The total amplitude Atot = VA^ -h in a (p, t)-plane for (P, 0 ) — (2, n7r/4) (n = 0,1, • • •, 7). 
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A. Initial amplitudes 


Let us see the initial amplitudes for the ingoing and outgoing waves. At the initial time t = 0, we have the ingoing 
and outgoing wave amplitudes 


Aq 


N + 12^aiarq^{p^ + q^)p 

D ’ 


B = 
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Aq IN — 128 a^a^g^(p^ + q‘^)p 

V D 


(34) 
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where 

N := (|a|^ + 6Aa‘fq‘^)p^ + 64g"^(|ap + )p^ + 64g^(|ap + a^)p^ + 64g^a^, 

D := (|n|^ + 64^^)^p^ + 256^'^(|(2|^ H- H- 64^^)p^ + 128^^(|(2|^ H- 4(2^Atot H“ 192^^)^"^ + 256^^((2^ H- 64^^)p^ + 4096^^^. 

As is shown in FIG El the disturbances for the total amplitude ^tot, which is proportional to the C-energy density, is 
localized in the neighborhood of the axis. 
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FIG. 2: Time-dependence of total amplitude Atot for {k,0,q) = (2,n7r/4,1) (n = 0,1,2,3): The blue-colored and red-colored 
graphs denote incident waves and refiectional waves, respectively, for t =, ±1, ±2, ±3. The violet-colored graphs in the case of 
n = 0, 2 mean that incident and refiectional waves entirely overlap. 


B. Asymptotic behaviors 

1. Timelike infinity 


Next we consider the asymptotic behaviors of the waves at late time t ^ oo. At t ^ oo, the metric behaves as 
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Here let us introduce the new coordinate 5 := 2 ; + so that djdz is also a translationally symmetric Killing vector, 
and then we can show that this asymptotic metric is that of the Minkowski spacetime. Therefore, both ingoing and 
outgoing waves vanish at late time. We have the asymptotic amplitudes 

Ac. Be. ^ + 0(^-3). (37) 

In particular, for = 0, the amplitudes behave in a different way as 

A c B c ^ + o(r^). (38) 

In this region, the polarization angles Ba and Ob for outgoing and ingoing waves behaves as 

tan6>A — tan6>B 1. (39) 

This shows that regardless of the values of the parameters, the x mode becomes dominate at late time. 


2. Spacelike infinity 

Let us consider gravitational waves near spacelike infinity p ^ 00 . In the limit of p ^ oc, we have the following 
asymptotic metric form 
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we immediately find that this spacetime is Minkowski spacetime with the deficit angle 
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The wave amplitudes behave as 
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For Ui 7 ^ 0, the polarization angles approach the values 

P - y|a|4 + IQajq^ 
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and for = 0, they behave as 
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3. Axis 


Furthermore, let us see the behaviors of the waves on the axis of symmetry p = 0. Near the axis, the metric behaves 


as 
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Using the coordinate z = z + Urcf), we have the metric 


ds^ - 




2\2 


4(^2 + q^Y + (tar - qaiY 


dz^ + 


4{t‘^ + q^Y + {t^r - q^iY 
4(g2 -(- t‘^Y 


p^d<Y + 


4(t2 




It is straightforward to show that the ratio of length of the circumference to the radial distance on the axis is 


lo^ y/^d(p 

hm % '' -— 

Jo V^dp 


= 27T. 
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This equation means that no deficit angle is present on the axis, which is in contrast to spacelike infinity. In the limit 
of p —> 0, the wave amplitudes behave as 
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The polarizations are 
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4- Null infinity 


At null infinity v ^ oo^ the wave amplitudes behaves as 
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where 


Ni = + 64afg®) - 128arai9Pi(9^+Pi)^ + 64(|ap + 0 ^) 5 "^^! + 64(|a|^ + a^) 5 ^pf + 64a^p}^, (53) 

Di = q^{\af + 64g^)^ + 256g^(|a|^ + aj. + Q4q'^)pf + 1024araiqp^{q‘^ - pj) + 384q'^(3|a|^ - 4a'J. + 64q'^)pl 

+256(a^ + 64g^)pf + 4096pf, (54) 

N 2 = \a\'^q‘^ + Q4ajq^ — 384araiplq{q^ + pI) + 192q‘^p‘J{\af — 5aj){q‘^ — Pi) + l^SOarttiq^Pi + Q4alpl‘^, (55) 
D 2 = q'^{\a\'^ + t)4q^)‘^ + 256g"^(|a|^ + aj 4- t)4q^)p\ + I924araiqp\{q^ - p\) + 384q'^(3|a|^ - 4a?r + 64q‘^)pl 

+256(a2 + 64^2)^12 (56) 


and Pi = J2u+ \/4u^ + q^. 


C. Faraday effect 


In this subsection, using the two-soliton solution, we study gravitational Faraday effect, which is a phenomenon 
that an outgoing (ingoing) wave amplitude corresponding to the + mode converts to an outgoing (ingoing) wave 
amplitude corresponding to the x mode due to the interaction with an ingoing (outgoing) wave with the x mode. 
Let us see how the + modes convert to the x modes while the corresponding gravitational waves are propagating 
along null rays from the axis p = 0 to null infinity v ^ oc. It is the most interesting to analyze it in the cases of 
— 16q^ — SaiQ > 0 and — 16^^^ + Saiq > 0 since both the polarization angles Oa and Ob completely vanish on the 
axis four times, i.e., at t = t±±, where t±± are defined by, respectively. 
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ar ± — 16g^ — 
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(57) 
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As is shown in FIGISI and Ob have time-dependence on the axis p = 0 and no x mode is present there just only 
at that moments. 

The upper-left graph in FIGH] shows that the pure + mode wave passing (t,p) = (t++,0) partially converts to the 
X mode wave and that its conversion comes to stop asymptotically ai v ^ oo. The ratio of the x mode to the -f mode 
monotonously approaches to a certain constant value. In this process, a complete conversion of + to x, or of x to -h 

does not occur. The upper-right graph in FIGH] illustrates that the pure + gravitational wave passing (t, p) = {t-\ _, 0) 

converts a little to the x mode and soon converts to the pure + mode. After that, a little conversion occurs and 
again becomes the pure -1- mode twice and Ob approaches to a non-zero value. Finally, its conversion comes to stop 
asymptotically at v ^ oc. The ratio of the x mode to the mode becomes a certain constant value. The lower-left 

graph in FIGH) shows the pure -j- gravitational wave passing (t,p) = {t _^,0) completely converts to the x mode and 

after that partially converts to the + mode. At 'T’ ^ cxd, the x mode becomes dominate. The lower-right graph in 

FIGH) illustrates that the pure H- gravitational wave passing (t, p) = (t_, 0) partially converts to the x mode. At 

V ^ oo^ Ob asymptotically approaches to a certain constant. 



FIG. 3: The behaviors of | sin2^A| and | sin2^B| for {ar,ai,q) = (10,0,1) on the axis 


Finally, let us consider the gravitational wave passing (t, p) = (tx, 0), where tx = aiqjar, when no H- mode is present 
at the axis. As shown in FIGEl fhe pure x mode partially converts to the + mode but x mode again increases, and 
at V ^ oo, the value of Ob becomes a constant. 


D. Time shift 

Let us investigate a time shift phenomenon as the non-linear effect, which means that a wave packet propagates 
at slower speed than light velocity. Basically following the analysis in Ref. M , where how to measure a time shift 
for gravitation solitons was proposed, we numerically analyze the asymptotic behavior of the wave packets at future 
null infinity v ^ oo and past null infinity u ^ —oo. In principle, we can find a time shift of the wave amplitudes 
by comparing its arrival time at future null infinity with that of a massless test particle starting off past null infinity 
at the same time. As is illustrated in FIGIHI an incoming massless particle propagating along = 0 from past null 
infinity u = —oo arrives at an axis p = 0 and after reflection, in turns, it propagates along u = 0 toward future null 
infinity v = oo. Let us see how slow a wave packet which has a peak near u = 0 and -i; = 0 is, compared with the 
massless particle. 

In FIG 111 fhe blue-colored graphs denote the amplitudes for ingoing waves near past null infinity hm^i^_oo 
and the red-colored graphs denote the amplitudes for outgoing waves near future null infinity lim^^oo B^/v for 
{k^O) = (2,n7r/4) (n = 0,1,2,3), where the amplitudes are multiplied by {—u)^ and because of apparent decay 
at null infinity. We can interpret these results as follows. Let us consider an incoming massless test particle starting 
from past null infinity and propagating along a null geodesic v = 0. The particle is reflected at an axis p = 0 and then 
it propagates to future null infinity along a null geodesic u = 0. An observer at past null infinity sees an ingoing wave 
packet earlier than an incoming radial photon, while at future null infinity he sees the outgoing wave packet after the 
outgoing photon. This shows a time shift phenomenon. Moreover, we would like to comment that for the very large 
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V-+ V_ 


FIG. 4: |sin2^A| and |sin2^B| along the outgoing null rays r^++, r^+-, u _foi* {ar^ai^q) = (10,0,1) 



FIG. 5: I sin2^A| and | sin2^s| along the outgoing null ray Ux for (ar,ai,^) = (10, 2,1) 


values, regardless of the difference between ^’s values, the wave packet always propagate slower than a massless test 
particle, as is seen in FIG. 0 It can neither collide nor split in the process. 


E. Collision, coalescence and split of solitons 

Besides the time shift phenomena, when k ^ \q\^ the ingoing and outgoing waves take various shapes depending on 
the phase, as is seen in FIG El As is seen in these graphs, at least, either of ingoing and outgoing waves can have two 
peaks. 

For n = 0 (^ = 0), there are two ingoing wave packets, one with a small peak and one with a large peak near past 
null infinity, and two outgoing wave packets, one with a small peak and one with a large peak near future null infinity. 
This obviously shows that two gravitational solitons collide, which occurs near the axis p — 0, and then the larger 
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FIG. 6: An orbit of the massless particle which is propagating in an incoming radial direction along f = 0 and after reflection 
at the axis p = 0 is propagating in an outgoing direction along = 0. 




FIG. 7: Time shift for {k,0,q) = (1000, n7r/4,1) (n = 0,1,2,3). The blue-colored and red-colored graphs show the incident 
waves at past null infinity and reflectional waves at future null infinity, respectively. 


one of two solitons overtakes the smaller one. For n = 6 {0 = it j2), conversely, the smaller one collide with the larger 
one and then overtakes. 
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For n = 3 (6> = 7r/4) , there are two ingoing solitons at past null infinity but a single outgoing soliton at near future 
null infinity. This shows that two solitons coalesce (in reflection at the axis). In contrast, as seen for n = 9 (^ = 37r/4), 
a single wave packet splits into two wave packets. Such phenomena do not happen for other solitons, such as solitons 
of KdV equation. 



FIG. 8: Amplitudes of the ingoing and outgoing waves for (k^O^q) — (2,n7r/12,1) (n = 0,1, •••,11). The blue-colored and 
red-colored graphs show the incident waves at past null infinity and refiectional waves at future null infinity, respectively. 


V. SUMMARY AND DISCUSSION 

In this paper, using the Pomeransky’s inverse scattering method for a cylindrically symmetric spacetime and 
starting from the Minkowski seed, we have obtained the two-soliton solution which has two complex conjugate poles 
to the vacuum Einstein equation with cylindrical symmetry. As the one-soliton solution with a real pole (in our 
previous work 0 ), it has been numerically shown that the two-soliton solution (presented in this work) describes 
a gravitational wave packet with two polarizations that comes from past null infinity, is reflected at the axis, and 
returns to future null infinity. The one-soliton solution in 0 describes a shock wave pulse with infinite amplitude 
propagating at light velocity, which yields null singularities, but the two-soliton solution is entirely free from such a 
singular behavior. This fact itself should not be surprising because in the previous work a two-soliton solution 
with complex conjugate poles and without any singularities has been constructed. 

In this work, using the two-soliton solution, we have studied non-linear effect of cylindrically symmetric gravitational 
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waves, focusing particularly on (i) gravitational Faraday effect, (ii) time shift phenomenon and (iii) collision process 
of two solitons: 

(i) The polarization angles Oa and Ob of gravitational waves on the axis have time-dependence. In particular, if 

— 16q^ — SaiQ > 0 or — 16g^ + Saiq > 0, at the times t = t±± the x mode completely vanishes on the axis p = 0 

and the + mode only is present there. In this case, we have studied how the pure + mode on the axis converts to the 
X mode while it is propagating along the null rays u = u±±. It is shown that in any cases, the polarization angles 
asymptotically a certain non-zero constant, which means that both modes are present at future null infinity. 

(ii) The time shift we say here is a phenomenon that a wave packet of a gravitational wave propagates at slower 
velocity than light. This is slightly different from the context used in the field of usual soliton theories, where this 
term (which is also called phase shift) is used to mean that when two solitonic waves collide, each position shifts as 
compared with when it propagates alone. In the case of cylindrical gravitational solitonic waves, it is evident that this 
phenomenon is due to self-interaction of a gravitational wave when an ingoing cylindrical wave reflects rather than to 
the interaction by a collision of two solitonic waves since in a region far from the axis, this gravitational soliton seems 
to propagate at light velocity. 

(iii) For the two-soliton solution in this paper, we have clarified that two gravitational solitons can coalesce into a 
single soliton and also a single soliton can split into two by non-linear effect of gravitation waves. Such phenomena 
cannot be seen for solitons of other integrable equations such as solitons of KdV equation. For the KdV equation, 
when two solitons traveling in the same direction collide (the amplitude is not simply addition of the two individual 
solitons), each one soon separates from each other and then asymptotically comes to approach the same shape of a 
wave pulse as before the collision. 
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